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Abstract. We show that a random graph studied by loffe and Levit is 
an example of an inhomogeneous random graph of the type studied by 
BoUobas, Janson and Riordan, which enables us to give a new, simple, 
proof of their result on a phase transition. 



1. Introduction 

loffe and Levit have recently introduced and studied a new random 
graph model called a quantum random graph, more precisely a quantum 
version of the classical Erdos-Renyi random graph, because of connections 
with a quantum model. (A random-cluster type modification of the graph 
studied in [1] and here is a Fortuin-Kasteleyn representation of a mean- 
field Curie- Weiss model in transverse magnetic field, see and [1]. The 
physical interpretations of the parameters /3 and A below then are inverse 
temperature and magnetic field strength.) 

The purpose of this paper is to show that this random graph can be 
regarded as an instance of the general inhomogeneous model studied by 
Bollobas, Janson and Riordan j^l, and that the general results on a phase 
transition in jH] enable us to give a new, simple, proof of the main result by 
loffe and Levit 

The random graph model defined by loffe and Levit ^4 has three param- 
eters, P G (0, oo), A G [0, oo) and n G N, and we will denote the resulting 
random graph by Qn^- (loffe and Levit consider also the case /3 = oo, 
which gives an infinite random graph, but we will not do so.) We will keep 
P and A fixed and let n — > oo. The parameter n then measures the size of 
the graph; more precisely, as will be seen below, the number of vertices is 
random but roughly proportional to n; note, however, that typically n is 
not exactly the number of vertices. 

To define Qn^, following j^, let S/s be a circle of length j3, i.e., the interval 
[0, P] with the endpoints identified, and consider n disjoint copies S^, . . . , §^ 
of Sp. Next, let there be a Poisson process of "holes" on each circle S^; the 
Poisson processes have intensities A and are independent. Let Mi be the 
number of holes on S^; thus Mi ~ Po(A/3). The holes split each circle 
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into open intervals if, k = 1, . . . , M*, where M* := max{Mj, 1}, and we let 

these intervals be the vertices of Qn^- (If Mi = 0, so there are no holes on 
S^, we regard the entire circle as an interval.) In other words, the vertex set 

of Qn'^ is {/f :l<i<n, l<k< M*}; the number of vertices in Qn^ is 

thus ELi^r- 

Finally, for each unordered pair (i, j), we consider a Poisson process Lij 
on S^, with intensity 1/n and independent of everything else; for each point 
X in Lij we consider the corresponding points in and and take the 
two intervals l!^ and ij containing these points; we add an edge between 

these two intervals (regarded as vertices in Qn^). Note that this yields a 
multigraph; there may be multiple edges. For the purpose of this paper, all 
multiple edges may be ignored; more precisely, we may merge each group of 
parallel edges into a single edge. Moreover, the expected number of multiple 
edges is easily seen to be 0(1), more precisely, it is less than /3^/4, so they 
may be ignored for almost all other purposes too. 

Note that when A = 0, the vertices are the n circles S^, and Qn^ is the 
classical random graph G{n,p) with p = 1 — exp(— /3/n) ~ /3/n. 

We state the main result by loffe and Levit [Jj in Section [2 Our for- 
mulation differs somewhat from |^, but the main features describing the 
phase transition of Qn^ are the same. We give our proof in Section [IJ after 
showing in Section |31 that Qn^ is an instance of the inhomogeneous random 
graph in ,3 . 

We denote the numbers of vertices and edges of a graph G by v{G) and 
e{G), respectively. Moreover, for a subgraph G of Qn'^, the vertices are 
intervals in some S^, and it is of interest to study the total length of these 
intervals, i.e., X^/gy(c') |/|; we denote this sum by i{G) and call it the length 
e{G) of G. 

All unspecified limits are as n ^ oo. 

Acknowledgement. I thank Geoffrey Grimmett and Oliver Riordan for 
interesting discussions. 



2. Results 

We first observe that the number of vertices in Qn^ is random and, by the 
construction above, given by v{Qn^) = Z^iLi '^hils the length £{Qn^) = 
j3n is deterministic. Moreover, the number of edges e{Qn^) ~ Po(/5(2)n) 
(if we include possible multiple edges). It follows immediately from the law 
of large numbers, cf. 1)3. 3() and (|3.5() below, that, as n — > cxd, 

7;(Q^'^)/n ^EM* = + (2.1) 
^(Q^'")/n = /3, (2.2) 
e{Q^^)/n p/2. (2.3) 
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Since the expected number of multiple edges is 0(1), (|2.;-i|) holds whether 
we count multiple edges or not. 

Denote the components of Qn^ by Ci, C2, . . . , in order of decreasing num- 
ber of vertices (breaking ties by any rule), with Cj{G) = if Qn^ has fewer 
than j components. The main result of loffe and Levit is that there is a 
phase transition similar to the well-known phase transition in the random 
graph G{n,p) |2], [H]: for certain values of (/3, A), all components are of or- 
der smaller than n, while for other values of the parameters, there is with 
high probability (exactly) one giant component with a positive fraction of 
all vertices. More precisely, we have the following result, essentially due to 
loffe and Levit but in our formulation and with some details added. 

Theorem 2.1. Let 

F(/3,A) := ^(1-e-^^) -/3e-^^. (2.4) 
A 

If F{(3,X) < 1, let vr := 0, while if F{(3,X) > 1, let ir > be the unique 
positive solution to 



(A + 7r)2^ ' A + vr 



Then 

y^Ci)/n ^ p := ^ (1 - e-(^+-)^) + e"^/^ (l - e^^j (2.6) 

£{Ci)/n ^ vr/?, (2.7) 
e(Ci)/n ^ C := vr(l - ^/2)/3, (2.8) 
while viC2)/n ^ 0, £{C2)/n ^ 0, e{C2)/n 0. 

For A = 0, we interpret ()2.4|) as the limiting value 2/? — /3 = /3; the results 
then simplify to standard results for G{n,p). 

Note that if F(/3, A) < 1, then p = ^ = ir = 0, and thus all components 
are small {op{n) vertices and edges), while if -F(/3, A) > 1, then there is with 
high probability one (and only one) large component with a positive fraction 
of the vertices, length and edges. 

The constant vr is by (|2.7|) and (|2.2j) the asymptotic fraction of the total 
length of Qn^ that belongs to Ci. Hence, vr equals the asymptotic probability 
that a given (or random) point on IJi belongs to Ci. This leads to the 
following corollary, due to loffe and Levit 1^. (Our tt is denoted p(/?, A) in 
Note that our formulation is somewhat weaker than theirs.) 

Corollary 2.2. The asymptotic probability that two random points on |J" 
(chosen independently and uniformly) belong to the same component equals 
TT^. The same holds conditioned on Qn^ too, with convergence in probability. 

To choose a random point on |J" means that we choose a vertex in Qn^ 
with probability proportional to its length. There is a similar result (with 
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limit p^) if we choose two vertices at random with the uniform distribution, 
see Theorem 3.14(i)]. 

Remark 2.3. In the proof below based on general results in jH], tt is calcu- 
lated from the survival probability in a multi-type branching process. Be- 
cause of the special structure in our case, vr can also be interpreted as the sur- 
vival probability in a standard, single-type, branching process, with offspring 
distribution mixed Poisson Po(r^) where = min(r, j3) and T ~ r(2, A^^), 
see 15 . The multi-type branching process describes asymptotically the pro- 
cess of finding new vertices in a component by successive exploration of the 
neighbours of the ones already found; the type is here the length of the 
vertex. The single-type branching process ignores the lengths; this can be 
done because when we reach a new vertex, on a circle not seen before, its 
length has the distribution F^, independent of what has happened before. 

To consider this single-type process simplifies some of the calculations 
below; for example, it is elementary to see that this process is supercritical 
if and only if F{P,X) > 1. We prefer, however, to use the setup in P] 
below and leave it to the reader to explore modifications involving also the 
single-type process. 

We can also obtain other results on Qn^ from for example the asymp- 
totic distribution of vertex degrees (a certain mixed Poisson distribution) 
and the typical and maximal distances between two vertices in the same 
component, where distances are measured in the graph theory sense of the 
minimal number of connecting edges. 

3. The general inhomogeneous model 

The general inhomogeneous random graph G^(n, k) is defined as follows; 
see 131 Section 2] for further details. We proceed in two steps, constructing 
first the vertices and then the edges. Note that n is a parameter measuring 
the size of the graph, and we are primarily interested in asymptotics as 
n — > oo. (In general, as in Section ^ n is not exactly the number of vertices, 
see (|3.3j) below. Moreover, there is no need for the parameter n to be integer 
valued in general, although it is so in our application.) 

A generalized vertex space V is a triple (5, /x, (x„)„>i), such that the 
following holds. 

(i) 5 is a separable metric space. 

(ii) /i is a (positive) Borel measure on S with < //(5) < oo. 

(iii) For each random sequence {xi,X2, . . . ,X]\f,J of A^^^ points 
of S (where Nn may be deterministic or random). 

(Formally, we should wr ^\ . . . , ) , say, as we assume no rela- 

tionship between the elements of x„ for different n, but we omit this extra 
index.) 

Let M{S) be the space of all (positive) finite Borel measures on S, and 
equip M{S) with the standard weak topology: i/„ — > iff j f dun j f du 
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for all bounded continuous functions f : S —>M.. Let 



n 
1=1 



where 6x is the Dirac measure at a; € 5; thus random element of 

M{S). We will further assume that 

(iv) fin — ^ fJ-, as elements of M(5). 
Recall that a set ^ C 5 is a fi-continuity set if A is (Borel) measurable and 



fJ-idA) = 0, where dA is the boundary of A. The convergence condition (iv) 
is equivalent to the condition that for every /^-continuity set A, 

finiA) := #{i <Nn:xie A}/n ^ //(A); (3.1) 

see Lemma A. 2]. 

Given a generalized vertex space V, we further assume that k is a sym- 
metric non-negative (Borel) measurable function on S x S. We let the 
vertices of the random graph G^{n, k) be the integers 1, . . . , Nn- For the 
edges we have two different versions. (See |3J Section 2] for further, differ- 
ent but asymptotically equivalent, versions). In both versions we assume 
that Xn is given and consider each pair of vertices with i < j sepa- 

rately. In the first version we add an edge between i and j with probability 
Pij := min|K(xj, Xj)/n, l}. In the second version, we instead create a multi- 
graph by adding a Po(K(xj, Xj)/n) number of edges between i and j. In 
both versions, this is done independently (given x„) for all pairs (i,j). 

In order to avoid pathologies, we finally assume 

(v) K is continuous a.e. on 5 x iS; 

(vi) K £ L^{S X S,fi X fj,); 
(vii) 

^ Ee{G^{n, k)) ^ \ I [ k{x, y) dfi{x) dfi{y). (3.2) 



As shown in [31 Remark 8.2], (vii) follows from the other assumptions if k 
is bounded and Var(A^„) = o(n^). 

Note that the number of vertices v{G^ {n, k)) is roughly proportional to 
n; more precisely, by (|3.1I) . 

= — = f^n[o) — > l^{S). (3.3) 

n n 

3.1. The lofFe and Levit random graph. We now show that the con- 
struction of Qn'^ in Section n is an instance of the general construction just 
given. For the space S we may choose the set of all open intervals in S^. 
(We regard the full circle §^ as an "interval", and thus an element of 5.) 
The first phase of the construction defines the vertices of Qn'^, with each 
vertex corresponding to an interval li € S. In the second phase, edges are 
added as in the second version above, with k(/, J) = |/n J| for two intervals 
/,JCS^. 
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The construction of edges a.s. does not distinguish between intervals of 
full length [3 obtained by from a circle with a single hole, and the full circle 
obtained when there is no hole at all. We thus identify the different intervals 
of length (3 and let the precise definition of 5 be 5 := 5i U 52, where 
Si := {S/3} is a singleton and ^2 is the set of all (open) intervals in of 
length less than (3. It is convenient to use the identification 

52 = {(x, ^) : X G S^, < £ < /?}, (3.4) 

where x is the left endpoint and ^ is the length of the interval. We can define 
5 as a metric space by taking the product topology on ^2 and letting both Si 
and ^2 be open in 5; we omit the standard construction of a corresponding 
metric. 

Remark 3.1. It might be more natural to define 5, as a topological and 
metric space, as the cone obtained as a quotient space of the product x 
(0,/5] by identifying all {x,j3) to a single point. (This gives a space that 
can be identified with an open disc of radius (3.) For our purposes, the 
construction is not very sensitive to the choice of topology on 5, and both 
these choices work well. 



Clearly, (i) holds. Moreover, if x„ is the collection of all intervals (in 



some order), then (iii) holds. 

Consider the Poisson process of holes on a single circle S/3. This yields 
a (finite) random family of intervals Ii, ... ,1m*-, i-e., a random subset of 
S. (In other terminology, a point process on 5.) We let ^ be the intensity 
measure of this random set; in other words, for every Borel set ^ C 5, ^{A) 
is the expected number of the created intervals that belong to A: 

^x{A) :=E#{i < M* : h G A]. 

In particular, since the number of intervals M* = max{M, 1} with M ~ 
Po(A/3), 

^i{S) = E M* = E M + P(M = 0) = A/3 + e"^^. (3.5) 



Clearly, (ii) holds. 

Since the construction of holes and intervals proceeds independently in 
the n circles S^, (|3.1|) follows, for every Borel set A'Z S, immediately by the 



law of large numbers for sums of i.i.d. random variables. Hence, (iv) holds 



As already remarked, the second stage in the construction of Qn'^ is as 
required, with k(/, J) = |/ n J[ for I, J £ S. Clearly, k is continuous and 
bounded on S, and thus the remaining conditions (v)[ (vi) and (vii) hold, 
using (2l Remark 8.2] for the latter. 

The random graph Qn^ is thus an instance of the inhomogeneous random 
graph in 0, and we may use the results of pj. 

Note that k{I, J) = for some pairs of intervals /, J, but k{I, S^) > for 
every I € S. Since /xdSg}) > 0, it follows that k is irreducible in the sense 
of Definition 2.10]. 
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To calculate the measure /u, note first that 

fi{{Si3}) = P(M* = 1) = P(M < 1) = (1 + A/?)e-^^; 

this is the point mass at Si. 

On ^2 we use the coordinates {x, i) as in (jSHJ. An interval I = {x,x + i), 
with < £ < /3, is obtained as a piece if and only if both x and x + i are 
points in the Poisson process 7i of holes, and further no point of Ti. falls 
in {x,x + i). The intensity of Tl is A, and given that there is a point at 
X, the gap L to the next point (if we unwind the cicrle Bp to a line) has 
an exponential distribution with distribution function P(L > £) = e~^^ and 
thus density function Ae~'^^. Consequently, the restriction of [i to 52 has 
the density A^e"'^^ dx M. (Another, but related, argument to see this uses, 
twice, the fact that the Palm process of a Poisson process equals the process 
itself.) 

We thus have the formula 

/" / d^ = /(S/3)(A/? + l)e-^^ + f f fix, dx dl. (3.6) 

Js Jo Jo 

Note that taking / = 1 yields dfi = e^^l^ + A/3 in accordance with ()3.5() . 

4. Proof of Theorem 12.11 

Since Qn^ thus is an instance of the general inhomogeneous graph G^{n, k), 
the results in connect the asymptotic behaviour of Qn to properties of 
the integral operator 

Tfix) = k(x, y)f{y) dn{y) 

on L'^{S,ij,). For example, one of the main results |3i Theorem 3.1(i)] says 
that there is a giant component in Qn'"^ if and only if ||r||j;^2(5^^-) > 1. (For 
simplicity, we write || ||x also for the operator norm of an operator on a 
Banach space X.) 

Lemma 4.1. The operator norm \\T\\i^2i^g^^-^ equals F{(3, A) defined in 1)2.4(1 . 

Proof. We exploit the symmetry of the construction under the group R := 
{rylygs^ of translations (or rotations, depending on your point of view) of 
the circle Sp, where we define Ty{x) := x + y (mod (3) for x,y G S^. Each 
Ty acts naturally on S (with G ^Si as a fixed point) and on L^(5,//), and 
it is obvious that T is invariant under R, i.e. that Try = TyT, y € S/3. 

Note that a function /(/) is i2- invariant if and only it depends on the 
length |/| only. Let L|j(5, /i) be the subspace of L^(<S, //) consisting of R- 
invariant functions, and note that the i?-invariance of T implies that T : 
L|j(5,/i) L|j(5, /i). Since k is bounded, T is a Hilbert-Schmidt operator 
and thus compact on Lp'{S, /i), and there exists a unique positive normalized 
eigenfunction ijj with eigenvalue ||r||j;^2(5^^). Since every translation Tj^(V') 
evidently is another such eigenfunction, Ty{ilj) = if) so if) ^ L\{S,^). (More 
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precisely, Ty{tp) = tp a.e. for every y, and it follows easily that = ip' a.e., 
where Ty{Tp') = tp' everywhere for each y.) Consequently, 

\\T\\ms,p.) = \\T\\Ll{s,^l)■ (4-1) 

Let S := {0,(3]. If g is any function on 5 = {0,(3], we let g denote the 
function / i— > 5(1^^1) on S. Thus g is /^-invariant, and every i?- invariant 
function on S equals g for some g. Further, let fi be the measure on S 
induced by /i and the mapping / i— s- |/| of S onto S; in other words, 

gdfi = j^gdfi (4.2) 
for any non- negative function g on S. Explicitly, by (|3.6|) . 

gdfL= [ g{\I\)df,{I) = g{(3){XP + l)e-^^+ l\{l)p\\-^' dt (4.3) 



Equivalents, djl = {\(3 + l)e-^^8p + (3X^6'^^ di. Clearly, the mapping 
g ^ g defines a natural isometry of Lj^{S, fi) and L'^{S, fi). 

If /, J € S, the average of K{I,Ty{J)) = |/ n Ty{J)\ over all translations 
Ty is |/| \ J\/P- Hence, if € L'^{S,(j), 

Tg{I) = ^ k{I, J)g{\J\) d/.( J) = ^ ^-^g{\J\) d/z( J) = ^ ^-^g{l) d^^)- 

(4.4) 

Define the kernel k on 5 by k(x, y) := xy/ (3 and let T be the corresponding 
integral operator. We then have, by (|4.4I) . Tg{I) = Tg{\I\), and thus 

T5 = (f ^^)^. (4.5) 

Hence, the restriction of T to L'j^{S,fi) is unitarily equivalent to T on 
L'^{S,jl), and thus, using (|4.1j) . 

Finally, the kernel k has rank 1, and Tg = (3~^{g,i)L, where i{x) = x. 
Hence, by ()4.3|) and an elementary integration, 





= /3(A/3 + l)e"^^ + / e'^X^e-^^di = F{P,X). 
Jo 



□ 



By Lemma [4.1l and Theorem 3.1], v{Ci)/n — ^ p for some /9 > 0, where 
p = if F{(3, A) < 1 and p > if F{f3, X) > 1. Moreover, by Theorem 
3.6], v{C2)/n — ^ 0; hence also i{C'i)ln because i{C2) < (3v{C'2), and 
e{C2)/n — ^ by ^ Proposition 8.11]. 
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To find the constant p, let, as in [2], Pk{x)-, for x E 5, be the survival 
probability for a multi-type Galton-Watson branching process with type 
space 5, starting with a single individual of type x and where an individual 
of type y has offspring according to a Poisson process on S with intensity 
K{y , z) (lp,{z) . Then, see Theorem 6.2], Pk{x) is the maximal solution of 

= l-e-^^''(") = l-exp(- / K{x,y)p^{y)y, (4.6) 

J s 

further, if ||r|| < 1 then Pk{x) = 0, while otherwise Pk{x) > for every 
X G <S, and there is no other non-zero solution. Then, 

p = j^Pnix)dpix). (4.7) 
Furthermore, by [SI Theorem 3.5], e(Ci)/n — ^ where 

Similarly, as an immediate consequence of [31 Theorem 9.10], 

£{C,)/n^ jjI\pM)dp{I). (4.9) 

It thus remains to find Pk{x), and to compute these integrals. 

By symmetry, Pk{I) depends on |/| only, and thus p^ = 4> for some func- 
tion (j) on S. Then 1)4.61) and 1)4. 5|) yield 

= 1 - exp(-r0) = (1 - exp(f (/.))^ 

and thus 

</> = 1 -exp(-f(/>); (4.10) 

in other words, (p is the corresponding function for the kernel k on S. Since k 
has rank 1, T(j) = ttl for some vr > 0, where as above l{x) = x. Consequently, 
(|1TTU)) shows that 

(/.(x) = 1 - e"^^. (4.11) 

Since T(l){x) = Jg ^(l){y) dfi{y), the relation Tcp = ttl, or T(f){x) = ttx, can 
be written 

vr = ^ |0(y) dA(y) = r ' ^ y(i - e""^) dA(y). (4.12) 

Evaluating this integral by (|4.3() . we find the equation (|2.5j) . The results 
quoted from Theorem 6.2] above imply that if -F(/3, A) < 1, then tt = 0, 
while if -F(/?, A) > 1, then vr > 0; moreover, there is no other positive solution 
of 1)2. 5|) . since otherwise we would have another positive solution of (|4.1U|) . 
Next, by (jiTj) and (g^l), 
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By (|4.11|) . (|4.H)) and an elementary integration, we find (|2.6|) . Further, by 
(lOl) and (jUnj, 

^(Ci)/n — !— 5- / x0(x) d/i(x) = /37r, 



which yields (ETJ. Finally, and (|0) yield 

c = ^ ^^^^ ~ ^^^'>^^y'>') "^^(^^ 

= ^ ^2 ^ d/i(x) ^ y dfi{y) - (^j^ x(l){x) d/i(x)^ 

which by (|4.12|) and the elementary J^xdjl(x) = j3 yields H2.8|) . 
Proof of Corollarv \2.iA The sought probability is, given Qn'^, 

P(same component | Q^^'^) = ^f/^f/ = E(^)'- (4-13) 
By 1)2. 7() . the term with j = 1 tends to vr^ in probability. Moreover, 

/3n / ~ Bn ^ Bn ~ Bn ^ 
i>i ^ i>i 

Hence, the conditional probability in H4.13() tends to vr^ in probability. The 
unconditional probability is the expectation of l|4.13|) . which converges to 
TT^ by bounded convergence. □ 

Remark 4.2. The proof shows that the largest component in Qn'^ asymp- 
totically behaves as for the random graph defined using the kernel k on the 
space iS with the measure fi. Nevertheless, these two inhomogeneous ran- 
dom graphs differ in other respects, for example in the number of triangles, 
see jni Section 17]. 
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